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We define (super) n -energy tensors for non-gravitational fields. The possibility of 
interchange of superenergy between gravitational and other fields is considered. 



1 The Bel-Robinson tensor and its properties. 

The Bel-Robinson (BR) tensofli is denned by 

q-af3Xp _ (jap\o Q0 /.( _|_ q apXa fi ^-j 

where C ap \ a is the Weyl tensor and * indicates the usual dual operation: 

C a /3\^ = (1/2) r] a p pa C p,y x ^ = (1/2) r\\p.paC a ^° ' . Its main properties are: 
A It is completely symmetric and traceless 

tj-apXp, rj-(cxf3Xfi) rj^OL Xp Q 

B For any timelike unit vector u, we can define the BR super-energy (s-e) 
density relative to u as Wt (u) = T a 0Xp,u a u^u x u il . For any u this is a positive 
quantity as follows from the expression Wt (u) = E pa E p<J +H prT H / " T > where 

E a x (u) = Capx^u^u^ and H a \ (u) =C a /3Xp u^u^ are-the well-known electric 
and magnetic parts of the Weyl tensor, respectivelyBuu. Furthermore, Wt (u) 
is zero for some u if and only if (iff) the Weyl and BR tensors vanish 

{3u such that W T (u) = 0} <=> T Q/3A ^ = <=> C afjXp = 0. 

C The scalar Wt (u) still satisfies a strongerj-iaequality, which has been called 
the dominant super-energy property (DSEP)En3: for all future-pointing vectors 
v, w, x, y, we have T a p\ p v a u)P x x y^ > 0. This is equivalent to the "dominance" 
of the completely timelike component of T in any orthonormal basis {e/r-j-, that 
is, Wt (e*o) = Toooo > |Xj/3A^| Va, (3, A, /i. This property is very usefulla for it 
can provide s-e estimates and bounds on the "growth" of the Weyl tensor. 
D Another fundamental property is that the BR tensor is "conserved" in vac- 
uum, in the sense that T is divergence-free 

R af 3 = Ag aP =^ V Q T"^ = 

where R a p is the Ricci tensor and g a p the metric tensor. 

"Dedicated to Professor Llm's Bel on the occasion of the Spanish Relativity Meeting 1998. 
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2 The Bel tensor: properties and decomposition 

The BR tensor is the s-e tensor par excellence, but it was immediately gener- 
alized by Bel to a tensor keeping most of the properties above and valid for 
any gravitational field. The Bel tensofl is defined by 

2B af3Xp = R apXa R^ a +^R* apXa *R*^ a +*R apXa *R^ a +R^ apXa R*^ a (2) 

where R a f3\^ is the Riemann tensor and .the duals act on the left or the right 
pair of indices according to their position]. The properties of B are: 
A B a(3Xp = B ()aXp = B a(ipX = B Xpaf} , B aXp - = 0, but B al3 a p ^ in general. 
B The Bel s-e density is defined by Wb (u) = B a " XfJ, u a Uf3U\u fJl , so that 
2W B (u) = (EE) 2 + (HH) 2 + (HE) 2 + (EH) 2 > where EE, EH, HE and 
HH are the electric-electric, electric-magnetic, magnetic-electric and magnetic- 
magnetic parts of the Riemann-..tpnsor, which arc 2-indcx symmetric tensors 
spatial relative to u defined bjQO (EE) a \ (u) = R a pXfiU^ u p , (HH) a \ (u) = 
*R* a f3\^u' 3 u tl , (HE) a x (u) = *R a px fl u f3 u p ' , (EH) aX (u) = R^apx^u^. I use 
the notation (A) 2 = A liV ... p A pv "^ > for any such spatial tensor. Again, Wb 
vanishes iff the whole Riemann and Bel tensors vanish too 

{3u such that W B (u) = 0} ^> B a(3Xp = R afJ x P = 0. 

C The DSEP also holds for the Bel tensor: -Boooo > \B a px^\ for all aj/LA^i 
in any orthonormal basis. The proof of this result is easier using spinoraffli3. 
D Finally, the Bel tensor has non-vanishing divergence in general 

where Jx^p = —J^xfj — — V M i?A/3- Thus, B is conserved if J = 0. 

The decomposition of the Riemann tensor into irreducible parts induces 
a canonical decomposition B af}Xp = T a/3Xfl + + M af3Xp where Q and 

M. are called the matter-gravity and the pure matter s-e tensors, respectively. 
SeeB for their definitions and properties. Here I only remark that Q is not 
typical (it does not satisfy any of the properties A-D above) and that A4 is a 
good s-e tensor for the matter content, it satisfies all the properties of the Bel 
tensor including DSEP, and the pure-matter s-e density W_m (u) (defined as 
usual) is positive and vanishes iff i? M „ = 0. Furthermore, the Bel s-e density 
decomposes as the simple sum: Wb — Wt + Wm for every u. 

It is noticeable that B has non-zero divergence in general but is divergence- 
free in the absence of matter (then B = T). Hence, the question arises of how 
to define s-e tensors for matter fields to see if the conservation of s-e can be 
restored (M is not satisfactory because it has no sense in Special Relativity). 
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3 Definition of (super) "-energy tensors for arbitrary physical fields 

Consider any m-covariant tensor t (Lll ...u m as an r-fold (ni, . . . ,n r )-form (with 
ni + . . . + n r — m) by separating the m indices into r blocks, each containing 
riA (A = 1, . . . , r) completely antisymmetric indices. This can always be done 
because, even if t A11 ... Mm has no antisymmetries, it can be seen as an m-fold 
(1,. . . ,l)-form. Several examples are: = is a simple (2)-form, while 
V ' pFfiv is a double (l,2)-form; the Riemann tensor is a double symmetrical 
(2,2)-form (the pairs can be interchanged) and the Ricci tensor is a double 
symmetrical (l,l)-form; a tensor such as t^p — t(^ y p) is a triple symmetrical 
(l,l,l)-form, etcetera. I shall denote £ Ml ... Mm schematically by £[ ni i . r nr ] where 
[ua] indicates the ^4-th block with ua antisymmetrical indices. Then, one can 
define the duals by using the * acting on each of these blocks, obtaining the 
tensors (obvious notation and 4 dimensions): 

t * * , t * * * 

[4— ni ],..., [n r ] [«i],...,[4— n r ] [4— m],[4— ii2],-,K] [4— ni],...,[4— n r ] 

There are 2 r tensors in this set (including t[ ni i . ...[n r ])- Contracting the first 
index of each block with u for all these tensors we get the "electric-magnetic" 
decomposition of t til ...u m - The electric-magnetic parts will be denoted by 

( EE ._E )[ ni -i],...,[n r -i], , C^^- H ')[ni],...,K_i-i][3-n r ] > etcetera. 

r T— 1 

These 2 r tensors are spatial relative to u and all of them determine t^ l flm 
uniquely and completely. Besides, t lll ___u im vanishes iff all its E-H parts do. Let 
us define the "semi-square" (t[ ni ],...,[n r l x t[n 1 ],...,[n r ]) by contracting all indices 
but one of each block in the product of t with itself 

- 1 

(t X t)Ai J ui...A r /n r = ~, _ t\ 1 p 2 ...p rll ,...,\ r a 2 ...u nr tp, 1 1 ...,Hr 

A=l {UA 1J ' 

The basic s-e tensor of t is dehncd as the sum of the 2 r semi-squares 
constructed with ir ni i ...r n i and all its duals. Explicitly: 

2T Al(tll ... Ar/ir {t} = (t[ ni ],...,[n r ] X t[n 1 ],...,lnr])x llil ...x rlir + + 

(t . , Xt , . ) (3) 

V [4-ni],...,[4-n r ] [4-ni],...,[4-n r ]/ M ^ 

A Expression ([}]) is a 2r-covariant tensor, symmetric on each (A^/x^-pair, 
and if t\ ni \ r n i is symmetric in the interchange of [riA]-blocks, then T{i} is 
symmetric in the interchange of the corresponding (A^M^O-pairs. Moreover, 
T{t} is traceless in any (A^^-pair coming from [2]-blocks (i.e. if ua = 2). 
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B The super-energy density W (u) = T\ lPl ...\ rPr {t}u Xl u Pl . . . u Xt u Pt is posi- 
tive for 2W (u) = {EE ...Ef + {E... EH) 2 + ... + (HH ...H) 2 >0, and 

{3u such that W(u) = 0} T Al/il ... VjUr {i} = <=> ^...^ = 0. (4) 

Actually, (it) is the sum of the squares \t tll ... Pm | 2 of all the components of t 
in any orthonormal basis which includes jL = Sq. 

C The DSEP holds for the general T{i}P which is one of the main justifica- 
tions and most important properties of definition (^|). 

Example 1: The massless scalar field . Let <j> be a scalar field satisfying 
V M V M = 0. First, by considering V M </> as the basic field, one can construct 
the tensor (^) and after expanding the duals we get 

1 

2 J 

which is the standard energy-momentum tensor. This tensor is identically 
divergence-free. Second, one can use the double symmetric (l,l)-form 
as the basic object, and construct the corresponding tensor (||) which becomes 
(this tensor was previously found by Bel in Special Relativity, see alsoliiJ) 

T a0 x p { VV0} = V Q V A 0V M + V^V^VaV^ - 

-ffa/3V A V^V M V P ^ - gx P V Q V p <pVf3 V p (f> + \^g a pgx P V^V^V 7 V p . (5) 

This is the basic s-e tensor of the scalar field. Its divergence can be easily 
computed and there appear several terms proportional to the Riemann and 
Ricci tensors. Therefore, it is a conserved tensor in flat spacetime. Now, one 
can go on and build the (super) 2 -energy tensor T q/ 3a^ti;{VVV</>} associated 
with the triple (l,l,l)-form V Q V^V^, and so on. This produces an infinite 
set of basic (super)™-energy tensors, one for each natural number n. The 
following fundamental result holds: "the basic (super)™-energy (tensor) of the 
scalar field vanishes iff the (n+l) th covariant derivative of <j> is zero". 

Example 2: The source-free electromagnetic field. Let F pv be a 2-form 

satisfying Maxwell's equations V p F v p = 0, V p F v p = 0. The tensor (|) for F is 

1 

which is the standard divergence-free energy-momentum tensor. To define the 
s-e tensor of the electromagnetic field one takes the double (l,2)-form V a F^„ 
as basic field and construct the corresponding expression (||) 

T q/3Am {VF} = V Q F Ap V /3 ^ + V a F w V p F x p - g aP V a F Xp V a Ff - 

-\g\ l X a F !jp V[ ) F'"> + \g aSi gx^rF a( ^ T F ap . (6) 



T^{F } = F Xp F/ - -gx P F pa F p ° 
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This tensor is not symmetric in the interchange of a[3 with A/i and is tracelesa 
in A/i. It does not coincide with previous s-e tensors for VF, such as that oftL3 
(which is simply T a p\^{X7F } + T\^ a p{V F}). Its divergence (with respect to 
the third index) is not zero in general, but it vanishes in flat spacetime. Again 
one can construct (super)™-energy tensors associated to the higher derivatives 
of F, and the important result is: "the basic (super) n -energy (tensor) of the 
electromagnetic field vanishes iff the n th covariant derivative of F is zero" . 

Example 3: The gravitational field. If the gravitational field is decribed 
by the Riemann tensor, then the s-e tensor (||) is exactly the Bel tensor (|J). 
Similarly, the corresponding expression (]|) for the Weyl tensor coincides with 
the BR tensor ([l]). One can also construct the s-e tensor for the Ricci tensor 
as basic field. This tensor has similar properties to those of the pure-matter 
s-e part M. of the Bel tensor, but it is not the same. In fact 

R 2 

2 Ma/sxti = T a \p^{R} + T afi fj\{R} - T a fjx^{R} + — (45 q (a# m )/3 ~ 9af)9\ii) 

where R indicates the trace-free Ricci tensor and R is the scalar curvature. 
However, it is remarkable that M. a /3\^ can be certainly obtained as basic s-e 
tensor for a field involving only the Ricci part of the curvature, namelyH 

MafiXfj, = TajSXfj, {R. .. — C....} 

where iL_ — C.... denotes the double symmetric (2,2)-form R a px^ — C a px^ 
(see e.g. □ for its explicit expression in terms of i? M „). This together with the 
Example 2 above makes it clear that the basic s-e tensor (||) for a given field 
is not the unique one with the same good properties. 

What is the arbitrariness in the definition of s-e tensors? To answer this, 
let us consider only the super-energy tensors (4-index s-e tensors). From the 
definition (||), the basic T a p Xtl satisfies T q/ 3 Aai = = T q/3aiA - The general 

s-e tensor formed with index permutations exclusively (to keep DSEP) is then 

which does not have any symmetry in general. However, it is automatically 
symmetric in the interchange of (af3) <-> (A/i) whenever the original T a px^ 
has this property, in which case we can redefine ci = c\ + c 6 , £2 = c 2 + c 5 , 
C3 = C3 + C4, £4 = C5 = £q = 0. In general, if the original T a px^ satisfies the 
DSEP, then the general T a /3\fj, will also satisfy the DSEP (and a fortiori the 
positivity of the super-energy and the essential property (^)) at least when the 
constants satisfy ci, C2, C3, C4, C5, cq > 0. Thus, subject to this condition there is 
a six-parameter family (reduced to a three-parameter one when the interchange 
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between pairs holds for Tapx^) of s-e tensors satisfying the fundamental DSEP. 
Furthermore, T a /3\fi is symmetric in a(3 iff c-i = C4 and C3 = C5 (or £2 = C3 in 
the special case), and symmetric in \[i iff C2 = C3 and C4 = C5 (respectively 
C2 = £3). This provides a three-parameter (resp. two-parameter) family of s-e 
tensors with the same symmetries as the basic one; Tapx^ is symmetric in the 
interchange of (a/3) <-> (X/i) iff ci = cq and C3 = C4 (resp., in general). In 
particular, this proves that there exists a two-parameter family of s-e tensors 
for the gravitational field satisfying the same properties of the Bel tensor, and 
another two-parameter family for the scalar field EJ. A possible way to avoid 
this arbitrariness is to take the completely symmetric part ^(apxn)- In the 
gravitational case, this provides the s-e tensor T( Q/3A/J ) = T a p\ tl + M^ a p\^ (mo 
matter-gravity part), which has been also put forward recently by RobinsoiE3. 

4 Interchange of super-energy between different physical fields 

The important question of whether or not super-energy has any physical reality 
or interest may find an answer by studying its interchange between fields (as 
energy does) and its possible conservation. One of the most striking features of 
the Bel tensor is that it is conserved in the absence of matter (when it coincides 
with the BR tensor). Similary, the s-e tensors (||) and (||) are conserved in the 
absence of gravity (in flat spacetime). The natural question arises: can these 
tensors be combined to produce a conserved quantity? 

To answer it, let us consider the propagation of discontinuities of the Rie- 
mann and other fields. The notation and conventions we use are those in 
RefE-3: E denotes a null hypersurface whose first fundamental form is g a b, 
(a, b, ... = 1, 2, 3 are indices in S). The normal one-form to E is (n^n M = 0), 
and a basis of tangent vectors to £ is denoted by e a (n^e^ = 0). Obviously 
n = n a e a and g a bn a = 0. As the null vector n is in fact tangent to E._Qp2 
needs to choose a vector field transversal to £ which is called the rigging oEj 
and denoted by I (n^ = 1). There are many different choices for £, but 
given any of them we can define wJJ as the three one-form fields completing 

with n M the basis dual to {£, e a }, that is, = 0, e^ui^ — 8%. We also put 
g ab = g^u^id^. Notice that this is not the inverse of g a b, which does not exist 
because E is null and the first fundamental form is degenerateEj. This last fact 
also implies that there is no canonical metric connection associated to g a b in S 
However, given any rigging one can define the so-called rigged connectiorolij 
in E by means of f £ c = w^e^V^e^, which is obviously torsion-free (T% c = T" & ). 
The covariant derivative associated to the rigged connection will be written as 
V. The discontinuity of any object v across E will be denoted by [v], as is 
customary. A very important point is that all the formulae appearing in what 
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follows are independent of the choice of the rigging vector I. 

The vacuum case. As is well-known, the discontinuity across £ of the Rie- 
mann tensor is determined by a symmetric tensor B ab on £ (in principle, six 
independent components), as follows 

[RapXft] = B ab (n a ajp - n/jw^) (n A w^ - n^w') . (7) 

However, in vacuum only two independent discpitinuities survive because 
g ab B ab = and n a B ab = 0. Then, one can shovkaO 

V a (B 2 n a n b n c n d ) = 0, 

where B 2 = B ab B ab = g ac g bd B ab B c d > 0. Moreover, the object inside brackets 
is directly related to the BR tensor by 

9 R 2 a b c d a b c d [srafiXfil 

Lo n n n n = w a w^WjU (I yl . 

Therefore, the discontinuity of the BR tensor is conserved along £ in vacuum. 

The Einstein- Maxwell case. If there is only electromagnetic field F M „ in the 
spacetime, and this is continuous, the discontinuities of the Riemann tensor 
take the same form (Q) and besides 

[VxF^] = nxf^, f^ = f a (%w° - n v uj^) , n a f a = 

so that there only appear±wo independent discontinuities in the first derivatives 
of F. It can be provecOEZl that 

V a {(B 2 + f 2 )n a n b n c n d }=0, (8) 

where f 2 = f a f a — g ab fafb > 0. The interesting thing is that this quantity 
can be related to the s-e tensor (^), if F^ propagates in vacuum, by 

2f 2 n a n b n c n d = w^w^ [T a ^{VF}] . 

Notice that in this case we have [B a f3\^] = [Tapx^], because the pure-matter 
part is continuous. Thus, if a Maxwell field propagates in vacuum, the sum of 
the discontinuities of the Bel tensor and of the s-e tensor for the electromag- 
netic field defined in M) is conserved along E. A version, of this fundamental 
result and of relation were obtained by Lichnerowic£3 back in 1960. 

The case with a massless scalar field. Assume finally that there only exists 
a massless scalar field (f> in the spacetime whose gradient is continuous so that 
the energy-momentum tensor and the pure-matter part of the Bel tensor are 
also continuous, [A4 a( 3A M ] = 0. Then, expression ([?]) still holds and furthermore 

[V M Vr/0] = V n^n v 
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which provides a unique discontinuity freedom. Then, one can arrive a' 



V a { (fciB 2 + k 2 V 2 ) n a n b n c n d } = 0, Vfci, k 2 

and where, analogously as before, if the scalar field vanishes at one side of E 

2V 2 n a n b n c n d = w^w^ [T Q/3A ^{VV0}] 

for the s-e tensor defined in (g). However, in this case, given that fci and k 2 
are arbitrary constants, the discontinuities are divergence-free separately. In 
other words, under the stated assumptions, the discontinuities of the Bel tensor 
and of the s-e tensor for the scalar field defined in ^) are conserved along £ 
independently. 

Similar studies can be carried out for the higher-order (super)™-energy 
tensors. The relevance and interpretation of all these results are under current 
investigation. 
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